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Abstract -We develop a generic method in Liouville space to describe the dissipative dynamics 
of coherent interacting quantum dots with adiabatic time dependence beyond linear response. 
We show how the adiabatic response can be related to effective quantities known from real-time 
renormalization group methods for stationary quantities. We propose the study of a delay time 
as a characteristic time scale. We apply the method to the interacting resonant level model 
and calculate the nonlinear adiabatic charge response to time-dependent gate voltages, tunneling 
barriers and Coulomb interaction. The dot charge delay time is found to be given by a unique 
expression in all cases, in contrast to the capacitance and the charge relaxation resistance. We 
discuss the observability of the effects in molecular systems and cold-atom setups. 



Introduction. — Adiabatic transport through quan- 
tum dots associated with a slow cyclic time dependence 
of the system parameters has generated a lot of inter- 
est in recent years, particularly in connection with quan- 
tum pumps [IHl] and mesoscopic capacitors [SJ0, see 
also Rcf. U. For nonintcracting systems the scattering 
formahsm is a powerful tool to describe the adiabatic 
response [Dll]. A central challenge in this field is the 
understanding of the influence of strong interactions as 
they typically occur in small quantum dots. Although 
general current formulas have been derived in terms of 
Green's functions [SUlOj. their evaluation is quite diffi- 
cult in the coherent regime at low temperature. Progress 
has been achieved in the perturbative regime of high tem- 
perature where it was shown that pure interaction 
effects can be revealed by the adiabatic response which 
would be covered by more dominant effects in the steady 
state [TT1[T2]. These studies also included the properties 
of the i?C-time in linear response [13]. In contrast, for 
interacting quantum dots at low temperature, a generic 
formalism providing the adiabatic time evolution in re- 
sponse to any parameter in linear or nonlinear response 
is not yet available. So far, quantum pumping has been 
studied for special models, like e.g. the 2-channel Kondo 
model in the strong coupling regime [9] , the Kondo model 
at the exactly solvable Toulouse point [2] , and the single- 



impurity Anderson model within slave-boson mean-field 
approximation [15] . In addition, the research on interac- 
tion effects in mesoscopic capacitors has concentrated on 
the special case of linear charge response to an external AC 
gate voltage by using the standard relation to the equilib- 
rium density-density correlation function. Here, the main 
object of interest was the charge relaxation resistance R 
defined by expanding the charge response 5Q in the exter- 
nal frequency w via 



5Q = {C + iujRC'^)5V 



(1) 



where Q{t) = SQe^^'^* denotes the charge, C is the static 
quantum capacitance and V{t) = (5Ve~*"* defines the ex- 
ternal AC gate voltage. For a single transport channel and 
provided that the Coulomb interaction is weak, a universal 
relaxation resistance R = h/2e^ was found in the coherent 
regime at zero temperature [5l[T6]. For interacting metal- 
lic dots and the single-impurity Anderson model the Shiba 
relation was shown to be a powerful tool to analyze R and 
its universality [TTHTO] . Using bosonization, the influence 
of Luttinger-liquid leads on R has also been studied [5D] 
and a numerical approach has been used away from the 
Fermi- liquid regime |21| . 

In this Letter we develop a general approach to deal with 
the adiabatic dissipative response, where the time scale of 



the external modulation i 



is much larger than the in- 
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verse of typical relaxation rates Fc, of a coherent quantum 
dot at low temperature including Coulomb interactions. 
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We show that the adiabatic response can be calculated 
very efficiently by using quantum field theoretical methods 
in Liouville space developed in Refs. [22j[23] and general- 
ized here for the case of time-dependent Hamiltonians (for 
approaches within Keldysh formalism see the review arti- 
cles [TIIH] and the recent development [23]). We provide 
a general relationship of the adiabatic response to effec- 
tive Liouvillians and vertices known from real-time renor- 
malization group (RTRG) in the stationary limit with in- 
stantaneous time parametrization, based on powerful tech- 
niques for the calculation of Laplace- variable derivatives, 
recently used within the E-Row scheme of RTRG [55] . As a 
consequence, our formalism is suitable to any model which 
can be treated by RTRG, which is applicable for many 
generic models with charge and spin fluctuations [H] . Re- 
cent applications of RTRG cover the Kondo model for 
both weak [2M5D] and strong [35] coupling, and the in- 
teracting resonant level model (IRLM) [3T]. Most impor- 
tantly, in contrast to previous research, our formalism al- 
lows for the adiabatic variation of any parameter in non- 
linear response, where no relation to equilibrium density- 
density correlation functions is possible and where certain 
identities like e.g. the Shiba relation are no longer appli- 
cable. Therefore, going beyond previous studies of the 
linear response to an external gate voltage, we also study 
the response to other parameters, like the tunneling cou- 
pling or the Coulomb interaction, which experimentally 
can be either realized intentionally, or indirectly induced 
by the gate voltage. We even cover the regime of non- 
linear response, motivated by recent works on mesoscopic 
capacitors in the nonlinear driving regime |32| . Instead 
of the linear response formula ([l} for the charge variation 
by an external gate voltage, we decompose the dynamics 
of any observable A, with a nonvanishing instantaneous 
contribution and its adiabatic correction in response to 
any parameter, as A{t) = A^^^{t) + A^°''^{t). The central 
quantity of our interest is the delay time scale ta for the 
expectation value A, defined by 

TA = \A^^^/A^'^\ , (2) 

which describes the delay of the full solution comparing 
to the instantaneous one. For A = Q in linear response to 
a time-dependent gate voltage, it is equivalent to the RC- 
time. In general ta can be quite different from typical 
relaxation times, depending on the observable, the type 
of excitation and its amplitude, and it is of fundamental 
interest to understand its dependence on interactions. 

We use our method to consider the IRLM with a 
single lead, which constitutes a minimal model for the 
mesoscopic capacitor with one single- particle level, where 
strong correlations are induced by a local Coulomb inter- 
action between the dot and the lead. Recently, the IRLM 
has been extensively used to study nonequilibrium trans- 
port through interacting quantum dots [3T ]l33ti36] , includ- 
ing the dynamics of the time evolution into the stationary 
state [31] . We calculate the nonlinear adiabatic charge re- 
sponse and the delay time tq, including renormalization 



effects of the tunneling enhanced by correlations. Impor- 
tantly, we find that the functional form of the charge delay 
time Tq is robust against the choice of the time-dependent 
parameter even in nonlinear response, whereas the capac- 
itance C or the relaxation resistance R get a more com- 
plex form when the tunneling or the Coulomb interaction 
are varied. Finally, we analyze further possible experi- 
mental implementations of the predicted results for the 
IRLM with time-dependent parameters, namely via the 
Anderson-Holstcin model in molecular electronics or via 
the spin-boson model in cold-atom setups. 

Method. — We start from a general Hamiltonian 
H{t) = i7,es + Ea tJ-a{t)Na. + i?dot(0 + V (t) of an inter- 
acting quantum dot coupled to noninteracting fermionic 
reservoirs with time-dependent chemical potentials fj.a(t) 
and a flat d.o.s. of width 2D via the coupling V{t). Gen- 
eralizing the Liouvillian approach of Ref. [22] to the case 
of time-dependent Hamiltonians, one flnds that the dissi- 
pative dynamics of the reduced density matrix p{t) of the 
dot can be described by the effective Liouvillian equation 



idtpit) = / L{t,t')p{t')dt', 



(3) 



where L{t,t') is the effective dot Liouvillian obtained by 
integrating out the reservoirs (h = 1). At the initial time 
to the total density matrix factorizes into an arbitrary dot 
and an equilibrium reservoir part. Since we are only in- 
terested in the asymptotic dynamics we set = —oo be- 
low. Following Ref. [12], the effective Liouvillian can be 
calculated diagrammatically, where each diagram of order 
0{V") consists of a product of vertices G{ti), ti > ■ ■ ■ > 
tji, with effective dot propagators Il{ti,ti-i) in between. 
In addition, the information of the Fermi distribution and 
the d.o.s. of the reservoirs is contained in time-independent 
reservoir contractions between the vertices. Using the for- 
mal definition G(t,t') = G(t)S(t - t' - 0+), we find that 
each term can be written in terms of a generalized convolu- 
tion in time space as {G oHo ■ ■ ■ oHo G){t,t'), where {A o 
B){t,t') = J^,dTA{t,T)B{T,t'). Introducing the partial 
Laplace transform A{t;E) = f_^dt'e'^^+'°^'^*-*'^A{t,t'), 
we get 



(Ai o ^2 o . . . o An){t; E) = e^Ai{t- E) . . . A^{t- E) = 



e'^.»'^-.^^^A,{h,E,)...Ar,{t^,E„ 



Ej—E.tu—t 



(4) 



The special differential operator V — id^ prescribes 
the energy derivative to act left to the time derivative. 
This rule is a natural generalization to Laplace space of 
analog identities in Fourier space used for gradient expan- 
sions in the Keldysh formalism [37] . Formally, it allows for 
the straightforward application of the Liouvillian approach 
to time-dependent Hamiltonians, with the difference that 
the exponential differential operator has to be taken be- 
forehand. In the adiabatic case, the exponential can be 
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expanded in dEdt ~ ^ ^ 1, leading to an expansion of 
the effective Liouvillian, 



L{t;E) = L'''\t;E) + L^''\t; E) + 



(5) 



Here, L^'^\t; E) denotes the instantaneous part, where the 
time t enters only parametrically via the external parame- 
ters, and L^°'\t; E) is the first adiabatic correction, which 
is linear in the time derivatives of the external parameters. 

Once the effective Liouvillian L{t\ E) is known up to the 
adiabatic correction, one can use it in eq. ([3|) which reads 



idtp{t) = {L o p){t; 0) = e^^-^f L(t; E)p{t) 



E=0 



(6) 



in the mixed (t; £')-representation. Expanding p{t) = 
p^^\t) + p'-'^\t) + . . . analogously to eq. ([5]), we find by 



comparing equal powers in the external frequency 



Tr = 1 



Trpf'') = 



(7) 

(8) 



In all arguments of L^*/"^ and dEL'^^\ E ^ has to be 
taken. From these equations the instantaneous density 
matrix p*^'^ (t) and the first adiabatic correction (t) can 
be determined. We emphasize that this approach is even 
applicable in nonlinear response in the amplitude of the 
external perturbations, i.e. only the time scale of the ex- 
ternal modulation needs to be large enough. Furthermore, 
it allows for an adiabatic modulation of any parameter of 
the Hamiltonian and is not restricted to a time-dependent 
gate voltage. 

The algebra of ([7]) and ([8]) can be easily evaluated for 
quantum dots with two accessible states. If additional 
conservation laws are present (as, e.g., charge conservation 
in the IRLM or spin-S^ conservation in the Kondo model), 
the nonvanishing matrix elements of the Liouvillian can be 
written as 



Eas.: 



'Lss,ss — i^s — ir/2 + isT' , Lgg^ss — , (9) 



where s = ± denotes the two states and s = —s. At E ~ 
we get from ([7]) and ([8]) that the instantaneous density 



matrix is diagonal pV' = Vr /T^'^ = 1/2 + sF't^'/Ff*) and 
the adiabatic correction fulfills 



(a) 



-p'^'' with 



r'(o 



(10) 



Below we use this result to evaluate the adiabatic response 
for the IRLM. 

Calculation of L'-'^K We now turn to the central 
issue of how to relate the adiabatic corr ection L'^''\t]E) 
to the instantaneous quantities known from RTRG. Us- 
ing (HI), we can formally write L^°-\t;E) = VL^'\t;E) = 
id^g^dl''^^'^L^'\t;E). The representation of L'^^{t]E) by 
its diagrammatic expansion specifies what "left" and 



"right" means for the derivatives with respect to E and t. 
As a first step, we represent the derivative OeL^'^^ by effec- 
tive vertices and propagators, using a method developed 
in Rcf. [35]. For generic models with two types of vertices, 
namely single (e.g. tunneling) and double (e.g. Coulomb 
interaction, exchange, etc.), we decompose it into two con- 
tributions in leading order and find 



dEL'^'\t-E) 



dEL^r\t-^E)+dELP{t;E) 



(11) 



The diagrammatic rules are explained in detail in Refs. [HI 
[25] . The single (double) circles represent full effec- 
tive single (double) vertices with effective propagators 
n'*'(t;_E) = E-L^'i{t-E) between (the Laplace variable 
is shifted by the frequencies and chemical potentials of all 
reservoir contractions crossing over the propagator). The 
left slash indicates Oe and the grey (green, color online) 
line represents the reservoir contraction given by the anti- 
symmetric part /(w) — ^ of the Fermi distribution function. 
All possible diagrams for c^^L*^*' can be classified by the 
number of lines over the propagator containing a deriva- 
tive. In the next step we perform the time derivative idt 
right to the energy derivative. The energy derivative is 
then shifted by partial integration to the reservoir con- 
traction (indicated by a (blue) cross) [25]. This yields 



L^r\t;E) 



L'u\t;E) 



^ y ^ + (I) A ij) 
O— — ^ + (!) A (!) 



o(5^^^i + o(!x>-»v-(5(!) 



(12) 
(13) 



where the right slash represents idt and the hat indicates 



the differential operator V = iS^^'c?"^'''. The frequency 
integral in both diagrams of (|12p is well-defined in the 
wide-band limit, so (|12p provides an explicit expression 
for the adiabatic correction containing the tunneling ver- 
tices in terms of renormalized vertices and propagators. 
In contrast, the frequency integrals in (|13p are logarithmi- 
cally divergent. We therefore take a second derivative with 
respect to E, yielding an RG equation for the adiabatic 
part, L^""* 
the double vertices 



jj {t]E), after partial integration. This contains 



dEL\f 



it,E)^\^^^-\^^^ . (14) 

Eqs. (|T^ and are the final results for the evalua- 
tion of adiabatic corrections of the Liouvillian in leading 
order, based on the instantaneous values of the renor- 
malized vertices and Liouvillian, which are obtained from 
RTRG0. For the adiabatic part of the propagator, appear- 
ing in the second diagram of ([T^ and ([Ti]) each, we insert 



^Provided that the frequency integrals converge, we note that 
our results can even be applied to a frequency-dependent d.o.s. in 
the leads. Otherwise, another derivative with respect to E may be 
required. 
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nW = n(')L('')nW + (ai5nW)(i(9tLW)nW. The first term 
does however not contribute to the adiabatic propagator 
in leading order. 

An interesting question is whether derivatives with re- 
spect to the Laplace and time variable commute in leading 
order, i.e. whether the adiabatic correction to the effective 
Liouvillian, eqs. ([T^ and ((HI), can be written as 



L 



(a) 



1 



it;E). 



(15) 



A similar relation was investigated so far only for noninter- 
acting systems |38j . To analyze its validity we introduce 
the complementing differential operator V = id'^^^^ d].'^^^ , 
where the energy derivative is taken right to the time 
derivative. Analogously to and ([T^ one finds in lead- 
ing order 



V'4\t;E) = - ^ i + d^.^^ , (16) 

dEV'4\t;E) = if^^- if^^ . (17) 

The inverted hat represents the differential operator V. 
Using idsdi = V + V , we can write V — ^^idsdi + 



i(2? — V) and, thus, the correction to eq. for Lp"-* 

(OeL^^^) is given by half the difference of ([T^ and 
(([13]) and (|17p ). We first address the second diagrams on 
the r.h.s. of these equations: their differences involve the 
expression 



-{V - v')ii^''> = n(') (-{V - v')L^A n(*) 



(18) 



+ i {(aBnW)(iatL(*))n« - n«(iati('))(9Bn«)} 

for the propagator. Here, the first term on the r.h.s. can 
be neglected in leading order, whereas the second one is 
only zero if the Liouvillian and its time and energy deriva- 
tive commute. For special cases this is indeed possible: it 
follows trivially for blocks where the Liouvillian is diago- 
nal, as e.g. the 2 x 2-block L'-^g = Sgs'^s^^ of eq. ([9|). 
For 2-level systems with conservation laws, sec eq. ([9]), it 

holds also for the block L),^ since the zero eigenvalue 
of the Liouvillian can be omitted in a propagator stand- 
ing left to a vertex averaged over the Keldysh indices [H] . 
Therefore, for this block one can replace the Liouvillian by 
its nonzero eigenvalue —iT^^\t; E) and the second term on 
the r.h.s. of ([T5|) is again zero. If this is given (or if the 
term can be neglected in leading order for certain models) , 
we can write the correction to eq. (jlSp generically as 



\idEdtLfit-E) 



dEL\;;>{t;E) = dE 



\idEdtL^^\t-E) 



(19) 



(20) 



From this result we observe another condition for the va- 
lidity of (|15p . namely that it should not matter whether 
the right or the left vertex is differentiated with respect to 
time, i.e. the two vertices should be equivalent. Whether 
this is the case, depends on the algebra of the model under 
consideration. For noninteracting systems one can take 
bare vertices and the reservoir indices of the two vertices 
have to be the same due to the reservoir contractions con- 
necting them. In this case the condition is fulfilled if the 
vertices do not depend on the level index of the dot states, 
e.g. through differently time-dependent coupling to differ- 
ent leads [33] ■ For interacting systems, the validity of 
is more restrictive. The renormalized vertices can be quite 
different from the bare ones and the vertices get an ad- 
ditional dependence on the Laplace variable E which is 
shifted by the chemical potentials of the reservoir lines 
crossing over the propagator standing left to that vertex. 
As a consequence, the two vertices are never equivalent in 
the presence of a bias voltage and correction terms def- 
initely occur for time-dependent voltages. As discussed 
below, for the particular case of the IRLM with one single 
reservoir, correction terms to eq. p3)) do not appear in 
leading order. 

Results. — We use the above developed method to 
analyze the response of a mesoscopic capacitor at zero 
temperature, described by the IRLM, where TJros = 
X]fe ^kO'\o'k describes a noninteracting reservoir with flat 
d.o.s. 1/ of band width 2D, Hdot{t) = e(t)c^c denotes 
a spinless single-level quantum dot with time-dependent 
level position e(i), and 



V{t) = 



Ml 

2tiv 

U{t) 



^ [c^ak + h.c.) 

k 

^ (etc -1/2) 4 



(21) 



fefe' 



is the dot-reservoir coupling with the bare time-dependent 
tunneling rate Yo{t) and the time-dependent dimensionless 
Coulomb interaction U{t). 

As shown above, we can evaluate the adiabatic response 
from eq. PU]) . where F'"^ and F'*^"' can be extracted from 
eqs. (Ull) and ^ together with the RTRG results for 
the instantaneous vertices and the Liouvillian derived in 
Ref. [3T]. For E = Q and leading order in [/, the results of 
Ref. ini read 



F = Fo 



D 



\e-iV/2\ 



2U 



BeT = i- 



C/F2 



2 + (E)2 



(22) 



F' = — — arctan — — , ScF' = — - 

TT F/2 ' ^ + 



2 , .r^' (23) 



where we have omitted the index (i) for the instantaneous 
quantities. Furthermore, the analysis in Ref. [21] shows 
that the Coulomb vertex is zero in leading order for the Li- 
ouvillian elements containing F and F'. Therefore eq. pop 
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is sufficient to evaluate F^^^ and r'*^"^. Inserting the alge- 
bra for the instantaneous tunneling vertices into eq. 
one finds that eq. ([TS]) is valid for the calculation of F''") , 
whereas for r'-"-' a correction term occurs proportional to 
dtU . This yields the total result 



2 e2 + (I)2 



2^ + 



(24) 



bmce a£;r,rw ~ o{u) we can neglect them in leading 
order in eq. ^ and, by inserting ([22]) to ([24]) into pO| . 
we find after a straightforward analysis for the charge re- 
sponse given by Q = ep+, 

g(^) = Com 4, Q(") = -R^dVdt^, (25) 
ei el 

where i?o = ^l^' ^^"-^ ~ f;^ e^+(r /2)^ • ^^'^ special 
case of linear response and when only e is varied with 
time, C = Co is the static capacitance and R ^ Rq the 
universal relaxation resistance, in agreement with ([T]). In 
contrast, when T is varied with intent or via an accidental 
(but experimentally unavoidable) gate voltage dependence 
of To or U, we obtain in linear response eq. ([T|) with 



C = Co 



where « f fa. 



e dV 



R = 



RoCo 



(26) 



9r9'^ In- . 



As a consequence. 



C and R are very sensitive to the variation of other param- 
eters, and logarithmic terms due to rcnormalization effects 
occur, if the Coulomb interaction U varies with time. 

In this general case, where also the renormalized F 
varies with time, we propose to analyze the time scale 
tq. From © and we get 



^0 



Q 



(a) 



r/2 



(F/2)2 



RoC'o 



(27) 



which is of 0(r^^) close to resonance e ^ F and of 0(F/e^) 
away from resonance. This result holds for any variation 
of e, Fo and U and is also valid in nonlinear response. In- 
teraction effects enter only weakly via the renormalized F 
given by ([22]) . Importantly, tq reveals the static capac- 
itance Co for a pure change of the gate voltage in linear 
response, with the advantage that tq can be determined 
in the presence of the variation of any parameter. 

The experimentally accessible time scale tq is thus an 
interesting quantity, which, for the case of the IRLM, is 
stable for the variation of any parameter in linear or non- 
linear response. We note that this time scale can vary 
quite drastically if other observables or other models are 
studied. E.g., the time scales tq and t/, when Q is re- 
placed by the current I = Q, are in general the same only 
in linear response. For the IRLM, the time scale tj shows 
similar logarithmic renormalizations in nonlinear response 
as they occur in C and R for time varying U. 



Realizations. — Several experimental realization of 
the IRLM exist, where the different parameters can be 
modulated in a controlled way. As we outline here, the 
applicability of the IRLM is not limited to the description 
of an interacting quantum dot, but allows the observation 
of the predicted effects for various physical systems. First, 
we show that the low-energy physics of the IRLM is equiv- 
alent to the one of the Anderson-Holstein model, as first 
predicted in Ref. [40] . This model is widely used in molec- 
ular electronics [41] and describes a single-level molecular 
quantum dot, having a vibrational degree of freedom with 
frequency il coupled linearly to the charge of the dot 



dot 



V 



eMC^'c + flb^b- Xn{b + b^)c''c, (28) 

(29) 



^^^(ctafc + /i.c.) 



The parameters ej\/, Fj\/, A and can be related to the 
effective parameters e, Fq and U of the IRLM. Applying 
a Lang-Firsov transformation |42] , the coupling to the vi- 
brational mode can be incorporated into the tunneling, 
leading to an effective level position, e = ea/ — A^fi, and 
tunneling rate, Fo = Fi^e"^ [42j|43]. If the vibration 
frequency CI is large compared to the other energy scales, 
the virtual intermediate states between the tunneling se- 
quences with one or more bosons can be integrated out. 
This produces terms with n ^ 2 lead operators in the ef- 
fective Hamiltonian. At large A, all cotunneling processes 
with n > 2 are exponentially suppressed, while the two- 
particle processes enter as an effective interaction. Hence, 
by integrating out all vibrational modes the Anderson- 
Holstein model can be mapped onto the IRLM with ef- 
fective interaction, U 
agreement with Ref. [1 



2^^, with n D. This is in 
where it was shown numerically 
that this formula has even a broader range of applicability. 
A modulation of the tunneling barriers is always accompa- 
nied by a modulation of the effective interaction U, since 
it is proportional to the tunneling rate F^/. In this case, 
our results predict that logarithmic corrections appear for 
C and R, whereas the time scale tq only depends on e 
and F via eqs. ([27]) and (|22)) . The Holstein coupling in the 
Hamiltonian allows for the observation of the dot charge 
via the displacement of the dot ~ (6 + b^). 

Finally, our results can be used to extract information 
on the relaxation behavior of systems described by the 
spin-boson model, namely, two-level dissipative systems 
connected to a large ensemble of oscillators 



E 



0Jgblb, + ^Y.9-iib, + bl)- (30) 



The spin-boson model can be implemented by a Bose con- 
densate of atoms trapped by a focused laser beam [44] . 
Such ultracold gases in optical lattices provide experimen- 
tal realizations for theoretical models with remarkably in- 
dependent tunability of parameters including the interac- 
tion strength, in contrast to usual semiconductor quan- 
tum dot setups. The system's behavior depends crucially 
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on the spectral coupling function. For the ohmic case, 
i.e. when the coupling constant obeys J2q9qH^ ^ ^q) = 
2aa;e-"/-°, the spin-boson model can be mapped onto the 
IRLM if a sa 1/2 (close to the Toulouse hmit) [45ll46|. 
The effective IRLM parameters are U = 1 — \f2a and 
To = A^/Z?. Changing the coupling of the Bose con- 
densate to the spin via a one generates a time-dependent 
effective interaction \J . The resulting response {S^ of the 
spin, identified with — p~)/2 in the effective IRLM, 
allows for the determination of the interesting time scale 
Ts^, given by ([27| . Especially in the biased case, where 
Ts^ ~ r/e^, this time scale is expected to differ signifi- 
cantly from typical relaxation rates F and F/2 for the di- 
agonal and nondiagonal components of the density matrix 
[311112] • 

Conclusions. — In this Letter we provide a generic re- 
lation of the adiabatic response to real-time RG results for 
the stationary case. The presented scheme allows for the 
variation of any parameter in linear or nonlinear response 
and provides criteria when the adiabatic correction to the 
Liouvillian can be calculated directly via energy and time 
derivatives of the instantaneous one. We suggest a delay 
time as an interesting time scale and show for the IRLM 
that its expression is robust against the choice of time- 
dependent parameters and their amplitude. We confirm 
the universality of the AC relaxation resistance, unless a 
time dependence of tunneling and interaction is present, 
revealing logarithmic renormalizations due to charge fluc- 
tuations. We proposed different setups in molecular elec- 
tronics and cold-atom systems to observe the effects ex- 
perimentally. 
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